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Abstract. The Mixmaster dynamics is revisited in a new light as revealing a 
series of transitions in the complex scale invariant scalar invariant of the Weyl 
curvature tensor best represented by the speciality index <S, which gives a 4- 
dimensional measure of the evolution of the spacetime independent of all the 
3-dimensional gauge-dependent variables except for the time used to parametrize 
it. Its graph versus time characterized by correlated isolated pulses in its real 
and imaginary parts corresponding to curvature wall collisions serves as a sort of 
electrocardiogram of the Mixmaster universe, with each such pulse pair arising 
from a single circuit or "complex pulse" around the origin in the complex plane. 
These pulses in the speciality index and their limiting points on the real axis seem 
to invariantly characterize some of the so called spike solutions in inhomogeneous 
cosmology and should play an important role as a gauge invariant lens through 
which to view current investigations of inhomogeneous Mixmaster dynamics. 



PACS number: 04.20.Cv 

1. Introduction 

The Bianchi type IX spatially homogeneous vacuum spacetime also known as the 
Mixmaster universe [1] has served as a theoretical playground for many ideas in 
general relativity, one of which is the question of the nature of the chaotic behavior 
exhibited in some solutions of the vacuum Einstein equations [2] and another is the 
question of whether or not one can interpret the spacetime as a closed gravitational 
wave [3, 4]. In particular, to describe the mathematical approach to an initial 
cosmological singularity, the exact Bianchi type IX dynamics leads to the BKL 
(Belinski-Khalatnikov-Lifshitz [5]) approximation involving the discrete BKL map 
which acts as the transition between phases of approximately Bianchi type I evolution. 
The parameters of this map are not so easily extracted from the numerical evolution of 
the metric variables [6, 7]. However, recently it has been realized that these parameters 
are directly related to transitions in the scale invariant part of the Weyl tensor [8, 9]. 
In fact this leads to a whole new interpretation of what the discrete BKL dynamics 
represents, by lifting its interpretation from transitions in apparently gauge-dependent 
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slicing geometry quantities to transitions in gauge-independent spacetime curvature 
invariants. 

For a given spacelike slicing of any spacetime, one can always introduce the scale 
invariant part of the extrinsic curvature when its trace is nonzero by dividing by that 
trace, thus reducing the information it carries to be equivalent to that contained in 
the scale invariant part of the shear tensor, which describes the anisotropy in the 
cosmological evolution of the spatial geometry of the time coordinate surfaces. In the 
expansion-normalized approach to spatially homogeneous dynamics, this corresponds 
to the expansion-normalized gravitational velocity variables [10]. This scale invariant 
extrinsic curvature tensor can be characterized by its eigenvalues, whose sum is 1 
by definition: these define three fimctions of the time parametrizing the foliation 
which generalize the Kasner indices of Bianchi type I vacuum spacetimes where they 
reduce to constants. We will refer to them as the generalized Kasner indices following 
Belinski and Francaviglia, who refer to the corresponding orthogonal spatial frame of 
eigenvectors as the generalized Kasner axes [11]. These quantities are defined for any 
timelike slicing of a completely general spacetime, where they will be functions of both 
the time and spatial coordinates. A phase of "velocity-tcrm-dominatcd" evolution 
[12, 13] (also referred to as Kasner evolution) is loosely defined as an interval of time 
during which the spatial curvature terms in the spacetime curvature are negligible 
compared to the extrinsic curvature terms. Under these conditions the vacuum 
Einstein equations can be approximated by ordinary differential equations in the time. 
These lead to a simple scaling of the eigenvectors of the extrinsic curvature during 
which the generalized Kasner indices remain approximately constant and simulate the 
Bianchi type I spatially flat spatially homogeneous Kasner evolution. 

The Weyl tensor can be also be repackaged as a second rank but complex spatial 
tensor with respect to the foliation [14] and its scale invariant part is determined by a 
single complex scalar function of its eigenvalues, a number of particular representations 
for which are useful. In particular the so called speciality index [15] is the natural 
choice for this variable which is independent of the permutations of the spatial 
axes used to order the eigenvalues, and so is a natural 4-dimensional tracker of 
the evolving gravitational field quotienting out all 3-dimensional gauge-dependent 
quantities [16]. In contrast with the eigenvalues of the expansion- normalized extrinsic 
curvature/shear which is a coordinate dependent measure of the deformation of the 
time coordinate surface, the eigenvalues of the Weyl tensor are spacetime invariants 
independent of the properties of the foliation and/or coordinates used to study a 
given spacetime. The speciality index is a particularly useful combination of those 
invariants which is automatically "expansion- normalized" and closely connected to the 
Petrov classification of the Weyl tensor, an important geometrical tool for classifying 
spacetime curvature which has found useful application first in numerically generated 
spacetimes where it is helpful in identifying gravitational radiation. 

During a phase of velocity-term-dominated evolution, the Weyl tensor is 
approximately determined by the extrinsic curvature alone, and hence the scale 
invariant part of the Weyl tensor is locked to the generalized Kasner indices exactly as 
in a Kasner spacetime, where the complex speciality index is confined to the interval 
[0, 1] of the real axis. Of course during transitions between velocity-term-dominated 
evolution where the spatial curvature terms are important, the generalized Kasner 
indices and the Weyl tensor evolve independently, but the transition between one set 
of generalized Kasner indices and the next is locked to a transition in the scale invariant 
part of the Weyl tensor. This idealized mapping, approximated by the discrete BKL 
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map between Kasner triplets, can be reinterpreted as a continuous transition either in 
the generalized Kasner indices associated with the expansion-normalized shear or in 
scale invariant part of the eigenvalues of the Weyl tensor. 

For spatially homogeneous vacuum spacetimes, the BKL transition is a 
consequence of a Bianchi type II phase of the dynamics which can be interpreted 
as a single bounce with a curvature wall in the Hamiltonian approach to the problem. 
One can in fact follow this transition in the Weyl tensor directly with an additional 
first order differential equation which is easily extracted from the Ncwman-Pcnrose 
equations expressed in a frame adapted both to the foliation and the Petrov type of the 
Weyl tensor. This type of Weyl transition in the spatially homogeneous Mixmaster 
dynamics can be followed approximately using the Bianchi type II approximation 
to a curvature bounce, leading to a temporally isolated pulse pair in the real and 
imaginary parts of the speciality index which represents a circuit ("simple loop 
around the origin") in the complex plane between the two real asymptotic Kasner 
points on the interval [0,1] of the real axis (a "complex pulse"), with profiles 
in time qualitatively similar to the pairs of spatial profiles present in expansion- 
normalized metric/shear conjugate variable pairs near the spikes recently found 
explicitly in spatially inhomogeneous Gowdy cosmologies by Lim [17], and numerically 
or qualitatively by previous authors [18, 19, 20, 21]. In fact the discontinuity 
in the asymptotic limit of the speciality index towards the singularity is a gauge 
invariant characterization of the presence of a "permanent true spike" in the example 
discussed in Appendix C. These particular spikes seem to describe inhomogeneous 
discontinuities that seem to arise in asymptotic spatially inhomogeneous BKL 
dynamics in neighboring curves entering a cosmological singularity. 

For the spatially homogeneous Mixmaster universe, the graph of the speciality 
index versus time serves as a sort of electrocardiogram (EKG) of the "heart" of the 
Mixmaster dynamics, stripping away all the gauge and frame dependent details of its 
evolution except for the choice of time parametrization, while the image of the complex 
curve which produces it is independent of the time as well. This single complex curve 
in fact represents the plot of the two real scale invariant scalar invariants against each 
other much like a velocity phase plane plot of velocity versus position during each 
single oscillation cycle of a mechanical system, here analogous to a single heart beat 
in this EKG analogy. Figs. 1-4 illustrate the relationship between the single complex 
pulse in the speciality index for Bianchi type II vacuum spacetimes and the temporally 
isolated pulses in its real and imaginary parts graphed versus time. Fig. 5 illustrates 
a sequence of three such pulses for a Bianchi type IX vacuum spacetime. The details 
of both the Bianchi II exact solutions and the Bianchi IX speciality index transitions 
are left to Appendix B. 

The decades of attention given to the Mixmaster dynamics focusing on the 
extrinsic curvature eigenvalue transitions have really masked transitions in the Weyl 
tensor. While all of the remarks made below for the spatially homogeneous Mixmaster 
spacetime extend to the class A Bianchi vacuum spacetimes and, apart from some extra 
considerations accounting for other sources of spatial derivative terms in the Einstein 
equations, probably also to inhomogeneous BKL dynamics [22, 17, 23], we will limit our 
detailed discussion to the Mixmaster spacetime alone. However, much of the interest 
in the Mixmaster spacetime has been its connection to inhomogeneous cosmological 
dynamics, so it is important to at least demonstrate how the speciality index fits into 
current investigations of spatially inhomogeneous exact and numerical cosmological 
solutions of the Einstein equations, adding a new perspective through which to view 




Figure 1. a) The pulses in the real and imaginary parts of the speciality index 
for a Bianchi type II Taub solution typical for Uao — » 1"^ where <Soo 1, here 
illustrated with Moo = 1-1, K = 1 so that ft; 0.9932, S-oc = 0.0597. 

b) The speciality index makes a large cardioid-like circuit in the complex 
plane starting and ending on the limiting points near and 1, close to the values 
for normal incidence scattering off the curvature wall which describes the locally 
rotational type D case. 

c) A closeup of the initial and final points in the complex plane of the previous 
figure, showing the increasingly horizontal approach to the real axis near the 
endpoint value 1 on the real Kasner interval. 



this work. In particular the speciality index reintroduces the spacetime point of view 
through the 4-dimensional Weyl tensor and its associated scalar invariants in contrast 
with the coordinate dependent space-plus-time 3-dimensional description primarily 
used in current approaches. In some sense this is more in tune with the real spirit of 
general relativity. 

We therefore briefly comment on a single example illustrating the spatially 
inhomogeneous case of G2-symmetric Gowdy spacetimes where exact solutions have 
been discovered by Lim [17] with a spike in the gauge-dependent metric variables, 
leaving most of the details to Appendix C. For these spacetimes, the speciality index 
is a function only of the two symmetry breaking coordinates r (Taub time) and a 
spatial inhomogeneity coordinate a;, and in the simplest case of permanent spikes which 
correspond to a single bounce off a single curvature wall (a "curvature transition") 
can be visualized as a vertical tubelike surface in 3-space with a vertical strip gap in 
the real limiting Kasner interval. This can be done by plotting the speciality index 
horizontally in the complex plane versus the vertical spatial inhomogeneity coordinate 
for all times, thus obtaining a "spatial pulse surface." For fixed x (horizontal cross- 
sections of this surface), the speciality index makes a circuit around the origin in 
the complex plane starting and ending on the real interval [0, 1] corresponding to 
the asymptotic Kasner indices characterizing the vclocity-term-dominated limits as 
r — > ibcxD, and its real and imaginary parts trace out a correlated temporally isolated 
pulse pair as functions of time, very similar to what occurs in a simple Bianchi type II 
spacetime (the model for a single bounce of a single spatially homogeneous curvature 
wall), a comparison made in Fig. 6. Instead for fixed Taub time r, the speciality 
index winds around the tubelike surface monotonically in x so that its real and 



EKG of the Mixmaster Universe 



5 




Figure 2. The pulses for a Bianchi type II Taub solution with Moo = ^/{u — l) ~ 
1.618034, if = 1 so that iSoo = 5— oo ~ 0.83275 and the crossing point of the 
previous cardioid-like circuit now touches the horizontal axis. 



imaginary parts instead realize a spatially isolated correlated pulse pair in that spatial 
variable. The "permanent true spikes" of Lim's exact solutions [17] may be seen as 
a singularity in this single unified spatial pulse surface describing the gauge-invariant 
and scale invariant spacetime information characterizing these solutions of the Einstein 
equations. Fig. 7 illustrates this pulse surface. 

The remaining "transient true spike" solutions involve two successive such 
curvature transitions, and hence a marriage of two such tubelike speciality surfaces 
glued together (two pulses). Fig. 8 shows a typical speciality plot for a fixed value of 
X for a transient spike case where the second pulse is just forming. For even more 
complicated G2 symmetric spacetimes exhibiting Mixmaster evolution, this speciality 
surface becomes much more complicated, with an infinite number of such pulses, with 
spikes arising in this process presumably in a more complicated fashion. 

2. Metric and curvature 

The Bianchi IX spacetime is a spatially homogeneous spacetime with isometry group 
SU{2) acting on the compact spatial hypersurfaces of the natural cosmological foliation 
by the cosmic time t [24]. The vacuum solutions of the Einstein equations for this 
spacetime lead to a natural orthonormal frame and dual frame which correspond to 
the diagonal form of the metric in a symmetry adapted frame 

ds2 = dt^ - [a{t) oj^f - [h{t) - [c(t) J'f , (2.1) 

where the three scale functions are functions only of the cosmic time i, and describe the 

time-dependent spatial anisotropy of the space sections. The signature H is used 

only so that the Newman-Penrose conventions of Chandrasekhar [25] can be followed 
below. The spatially homogeneous 1-forms i — 1,2,3, the natural orthonormal 
frame {cq}, a = 0,1,2,3 and associated Newman-Penrose frame j-^, n,m}, and the 
field equations are reviewed in Appendix A. 

In order to explore how the gravitational field evolves near the big bang or big 
crunch singularities ahc — > it is convenient to introduce a new time coordinate t 
called Taub time [26] by setting 

di = a6cdr (2.2) 
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Figure 3. Same as the previous diagram but with Uao = 2, = 1 so that 
u—oo = 1 and hence <Soo ~ 0.708, cS-oo = 1. Now the right horizontal intercept 
of the circuit has moved right to be exactly 1 and the tangent line to the circuit 
there is horizontal. 

which pushes these singularities for which abc — s- out to r — oo (big bang) or 
T oo (big crunch), assuming r increases towards the future, and hence is more 
useful than the cosmic proper time which compresses the interesting behavior near 
these limits. Let a dot denote the cosmic time derivative and a prime denote the Taub 
time derivative, related to each other by /' = abcf. 

The field equations are also more conveniently expressed in logarithmic metric 
variables a, (3 and 7, defined by a = e", b = e^, c ^ e'', so that 

a" - i [(e^^ - e^-'f - e*"] = , /3" - i [(e^^ - e^")^ - e^^] = , 

7"-^[(e^"-e^'3)^-e47] =0, (2.3) 
with the constraint equation 

a'/3' + /3'7' + I'a' + l(e2("+/3) + e^iP+f) + ^^(1+^^) „ i(e4" + + e^'') = . (2.4) 

These equations should be evolved backwards in Taub time if r increases towards the 
future, in order to study the BKL dynamics which aims towards an initial singularity 
customarily set att = 0,T— s-— 00. 

The matrix of mixed spatial extrinsic curvature tensor components with respect 
to the spatially homogeneous frame {K'^j) = — diag{a, is diagonal (like the 
metric itself) and one can define the time-dependent generalized Kasner indices as 
the diagonal values of the ratio of its components with its trace 

{K',)/K''k = diag(pi,p2,P3) = diag(a',/3',7')/(a' + /3' + 7') , (2.5) 

from which it is clear that the generalized Kasner indices represent the scale invariant 
part of this tensor's eigenvalues and satisfy pi + P2 + Ps = 1- This quotient 
extrinsic curvature tensor corresponds to the expansion-normalized shear variables in 
the EUis-MacCallum-Wainwright approach to the dynamics of spatially homogeneous 
cosmologies [10], related to the generalized Kasner indices by 

(E*,) = {K', - 6',K\/3)/{K\/3) = diag(3pi - 1, 3p2 - 1, 3p3 - 1) • (2.6) 
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Figure 4. The Bianchi type II pulses typical for u > 2 illustrated here for 
Moo = 2.73, X = 1 so that 5tx> ~ 0.500, S-cx> « 0.8003. 



The quadratic Kasner index constraint comes from the scale invariant 
Hamihonian constraint 

' ^K^^Y -pI+pI+pI-^ (2.7) 

^ (g2(a+/3) g2(/3+7) + g2(7+a)) _ 1 (g4a ^ ^4/3 ^ ^4^) 

(a' +13' + 7')2 ■ 

Under conditions where the right hand side is very smaU, the second Kasner constraint 
is satisfied: 

pI+pI+pI = 1- (2.8) 

In this case one can simphfy the the following expressions in the expansion-normalized 
shear tensor to find 

S^.E^fcE'^, = (3pi - If + (3p2 - 1)' + (3p3 - 1)^ = 81piP2P3 + 6 = 3det(S*,) , 
det(E^,) =(3pi-l)(3p2-l)(3p3-l) = 27piP2P3 + 2. (2.9) 



3. The Weyl tensor 

Using the unit normal eg to the foliation, one can introduce the complex (spatial) 
tensor (see Eqs. (3.52), (3.53), (3.62) of [14], note the misprint Qab instead of Qac) 

— QaP — {CaipS + i ~CQT,/35)egeg = Eap + iBap , (3.1) 

associated with the self-dual part of the Weyl tensor, where the spatial tensors Eap 
and Bap respectively denote the electric and magnetic parts of the Weyl tensor with 
respect to eg. Note that the sign before i in the expression for the dual of the Weyl 
tensor depends on the convention chosen for the oriented unit volume 4-form 77, chosen 
here following [14]: 775 jgj = 1- 

The spatial tensor Qap is symmetric and trace-free, i.e., satisfies the conditions 

Q"c. = 0, Qap^Qpc, Qo^pu^^O (3.2) 

and its nonzero components when expressed in any orthonormal frame adapted to 
eg form a complex symmetric tracefree 3x3 matrix (Qij), = 1,2,3). The 
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mixed form of this spatial tensor represents a complex linear transformation. The 
algebraic type of the matrix Q — (Q^j) provides an invariant characterization of the 
gravitational field at a fixed spacetime point. This Petrov classification based on the 
eigenvalues Ai, A2, A3 of Q, their multiplicities, and the number of linearly independent 
eigenvectors leads to certain canonical forms for the matrix Q called the normal forms 
which are listed in [14] . At most two eigenvalues are independent due to the tracefree 
condition on the matrix Tr Q = Ai + A2 + A3 = 0. This classification is insensitive to 
the scale of the eigenvalues, and so only depends on their scale invariant part. 

The most general Petrov type I corresponds to Q being diagonalizable, and a 
canonical or normal orthonormal frame for this type is one in which the matrix is 
diagonal. This is the case for the natural orthonormal frame introduced above for the 
Bianchi type IX metric (see Appendix A), where the type I normal frame conditions 
tpi = ''Pa = i'o — V'4 lead to the following expression from Table 4.2 of [14] 

Q = diag(Ai, A2, A3) = diag(V'2 - ipo, V'2 + V'o, -2V'2) • (3.3) 

The Petrov type D occurs as the subcase Ai = A2 and permutations, conditions which 
if maintained in time correspond to the Bianchi type IX Taub spacetime. The scale 
invariant part of the Weyl tensor can be expressed various ways. For example, one 
can choose some ratio of the two independent Weyl scalars like z = 'ipi/tpQ ov one can 
pick one of six possible permutations of the eigenvalue ratios like = /U23 = A2/A3 
and get 

Q = A3 diag(-(l + n), IX, 1) = V'o diag(z 
which leads to the relations 

_ A2 _ V'2 - V'O _ 1 + Z~'^ _ — ^1 + ^2 

^ A3 ■tl>2 + V'O 2 ' ^ V'O A2 - Ai 

Note that the choice of z rather than some other ratio constructed from the three 
eigenvalues is also arbitrary, and subject to a similar action of the permutation group 
like fj, [8, 9]. 



l,z + l,-2z), (3.4) 



1 + 2^ 



(3.5) 



4. The speciality index 

The Petrov classification distinguishes between algc;braically general spacetimes (type 
I) and algebraically special ones (types D, II, N, III and O) depending on the 
degeneracy of eigenvalues and eigenvectors associated with the complex matrix Q. 
One way of determining whether a given spacetime is general or special involves 
evaluating the scalar invariants determined by the trace of the square and cube of 

Q 

7 = ^TVQ2 = l(Af + Xl + \l) = -Lc«^,C^*„^ , (4.1) 

J = ^TrQ3 = i(A? + A^ + Ai) = iAiA2A3 = j^C^'^sC''' ^.C'^'' . 

For Petrov types I and then D and II in the top two levels of the Petrov hierarchy, 
both invariants I and J are nonzero, while both vanish for types O, N and III in the 
lowest level of the three level pyramid reprsenting the Penrose specialization diagram 
(Fig. 4.1 of [14]). Any algebraically special spacetime satisfies = 27 , so except 
for most degenerate types O, N and III in the Petrov hierarchy, one can introduce 
the so called "speciality index" [14] which is a well defined function of the eigenvalues 
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Figure 5. a) The pulses in the real and imaginary parts of the speciality index 
plotted for several consecutive bounces in a Bianchi type IX spacetime. Each 
collision with a curvature wall generates a pulse in the real and imaginary parts 
of the speciality index, corresponding to a circuit in the complex plane beginning 
and ending on the unit Kasner interval of the real axis. 

b) The speciality index plotted as a sampled point plot rather than a smooth 
curve for these same three consecutive bounces. 



''^11-^2,^3 by taking the dimensionless ratio (which is therefore independent of the 
overall scale of the eigenvalues, i.e., scale invariant) 



S 



27J^ 



6 



(A3 + A3 + A3)2 _ 9 (AiA2A3)^ 



/3 (A? + Ai + Ai)3 4(Af + Ai + Ai)3 



(4.2) 



having the value 1 for the algebraically special types D and II and which is obviously 
invariant under any permutation of the eigenvalues. Even though it may not be well- 
defined for Petrov types O, N and III, a well-defined limit may exist within a family 
of special spacetimes, as occurs in the Bianchi type I case of the Kasner family of 
vacuum spacetimes. 

Re-expressing the speciality index in terms of either of the two scale invariant 
variables z or /i one finds 

27 fJ-Hl+txf Z2(z2_l)2 



S 



(4.3) 



4 (l+/, + /,2)3 (^2 + 1/3)3 • 

This object is invariant under the permutation group and hence does not depend on 
the particular ordering of the spatial axes like any of the possible choices for fj, or z 



EKG of the Mixmaster Universe 



10 



do. It is the only independent scale invariant combination of the scalar invariants of 

the Wcyl tensor which is also invariant under a permutation of the eigenvalues. 

For Bianchi type I dynamics where the generalized Kasner indices are constant 
and satisfy the additional quadratic identity, there is only one independent parameter 
which can be taken to be the scale invariant ratio u = U32 — P3/P2, in terms of which 
one has the Lifshitz-Khalatnikov Kasner parametrization [27] 

-u 1+u u(l+u) ,^ 

1 + U + U-' 1 + U + U-' 1 + u + u^ 

One finds that /z = u is then real [16, 8, 9], so that z and S are also real, and the 
speciality index has zero imaginary part while the real part is the same expression in 
u as in /i. Indeed the speciality index is then explicitly 

27 27 u'{l + uf 27{w' - if 
S = Sk = — -rPlP2P3 = -rj:^-, ; — = / o 1 q\3 ' ('^•^) 

where u> = 2'u + 1 is the real value of — = 2/i + 1. a Kasner parameter found 
useful by Lim [17]. Using Eq. (2.9), this leads to the following relations with the 
expansion-normalized shear tensor 

det(S%) = 2 - ASk , S'jSJfeS'^i = 6 - 12Sk ■ (4.6) 

One easily sees that in this case the speciality index is confined to the real interval 
< Sk < 1, equaling 1 only when two of the three Kasner indices coincide and are 
nonzero (type D) and equaling only when two of the three Kasner indices coincide 
and are (type O) . As discussed in [9] , the symmetric group of permutations of triplets 
of numbers acts both on the eigenvalues of the extrinsic curvature and of the Weyl 
curvature, and hence also on the speciality index, which is invariant under this action. 
This freedom is useful in adapting the variables /x or 2: or u to transitions between 
velocity-term-dominated phases of the evolution. Note that the evolution equation 
(A. 7) is easily converted into one for jj, or any of the possible choices for the scale 
invariant Weyl scalar invariant variable, in terms of which the speciality index can 
be expressed if this differential equation is integrated numerically in parallel with the 
usual variables. 

The speciality index is not only a scale invariant scalar invariant of the Weyl 
tensor, but during velocity-term-dominated phases of the evolution it also determines 
the scale of the scalar invariants of the Weyl tensor in a simple way. The simplest 
combination of the three Weyl eigenvalues which is invariant under the permutation 
group and is representative of their scale is their product T = A1A2A3. For the 
standard Bianchi type I metric variables where constants have been absorbed into the 
definition of the spatial coordinates 

(a, b, c) = (if 1 , ,tP') = (ef ef=^, e^^^) , abc = fPi+f^+fs = t = e'' , (4.7) 

this scalar is simple [28] 

(Ai,A2,A3) =t-Vp2P3 (—,—,— ) =t~^-^s] (—,—,— 

\Pl P2 P3/ V 27 7 \pi P2 P3 

T = A1A2A3 = {pw2Pzft-^ = (-^^5^ • (4-8) 

It scales by a power of the square root of the absolute value of the spatial metric 
determinant abc = t, with a proportionality constant determined by the speciality 
index itself. Note that all of the algebraic scalar invariants of the Weyl tensor are 
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a) b) 

Figure 6. a) A plot of the speciality index for constant values of x ^ for one of 
Lim's Gowdy true spike solutions, shown here for the parameter values w = 0.1, 
Qo = 1, A2 = 1. The cardioid-like circuits lying mostly in the upper half plane 
correspond to the positive values x = 10^ for j = —5, —4, —3, . . . , 2 from inside to 
outside representing selected spatial points on a logarithmic scale, revealing how 
the circuit expands with increasing a; > 0. The corresponding curves mostly in 
the lower half plane are shown for the sign-reversed values of x and arise by a 
simple reflection across the real axis under the sign change. All circuits for x ^ 
start at the common left Kasner point and end at the common right Kasner point, 
but the single curve for x = degenerates to line segment from the left Kasner 
point to an isolated Kasner point further to the left on the unit interval [0, 1]. 

b) For comparison, the circuit (lower cardioid-like curve, compare with 
a; < curves to left) and its complex conjugate circuit (upper cardiod-like curve, 
compare with x > curves to left) to correspond to time reversal for x > for the 
corresponding Bianchi type II Taub solution seed plot with the same asymptotic 
Kasner states Moo = 0.55, tioo — 1 = u = —0.45 , where w = 2u + 1. 



determined by these two complex numbers S and T, although in the larger context 
of Kasner transitions in Mixmaster dynamics one must also take into account the 
contribution from the change in the rescaling of the spatial axes to achieve the Kasner 
form of the scale factors. 

In the Mixmaster context, the Kasner expression Sk as a function of the 
generalized Kasner indices or their related Lifshitz-Khalatnikov parameter extends 
to a time-dependent function on the spacetime locked to the trace of the cube of the 
expansion-normalized shear tensor. It agrees with the speciality index appromimately 
only during velocity-term-dominated phases of the evolution. This shear quantity 
expressed in terms of the variable u has been used by Garfinkle [21] in the larger 
context of Gowdy spacetimes. 

To summarize, the speciality index is a simple scale invariant scalar invariant 
of the Weyl curvature tensor which, unlike any of the possible /i or z variables 
whose values are shuffled around six different regions of the complex plane by the 
permutation group [9] , is independent of the ordering of the spatial axes and has real 
values confined to the closed interval [0, 1] during velocity-term-doniinated phases of 
the evolution. During a transition between velocity-term-dominated phases of the 
evolution in the BKL dynamics induced by a collision with a spatial curvature wall, 
it exhibits a complex pulse which joins its initial and final real values. In contrast, 
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a) b) 

Figure 7. Complex plots of the speciality index <S (horizontal) a) versus the 
spatial coordinate x > in the vertical direction and b) versus the Taub time t 
in the vertical direction for the Lim's Gowdy true spike solution of the previous 
figure. 

a) Ten horizontal cross-sections x = 1 ... 10 are shown, together with the 
two common asymptotic Kasner values (the two vertical lines through the real 
axis) which are the beginning and end points for any isotemporal curve, like the 
typical isotemporal curve r = 3 shown in the figure, which winds around the 
tubelike surface with a slowly expanding cross-section as x increases. For a; < 
the symmetry S{—x) = S{x) reflects the horizontal cross-sections across the real 
axis, both of which degenerate to a real line segment for a; = 0. 

b) From top to bottom the curves x = 1, 2, 3, 4, 5, 10, 100 are shown, together 
with the vertical asymptotic Kasner lines. Each experiences a pulse in a limited 
time interval which is revealed as spiky behavior in the real and imaginary parts 
of cS. For a; < the symmetry S{—x) = S{x) reflects these curves across the real 
axis. 



all of the variables used to describe the evolution itself are functions of the frame- 
dependent 3-diinensional metric or extrinsic curvature tensor components. In other 
words, the speciality index is in many senses a privileged spacetime curvature scalar 
which tracks the evolution, and whose graph (real and imaginary parts) provides a sort 
of electrocardiogram of the heart of the Mixmaster universe independent of the body 
of 3-dimcnsional gauge-dependent variables which house it, figuratively speaking. 

Note that we can also establish a correspondence between the speciality index of 
the exact Bianchi type IX evolution and the approximate one defined in terms of the 
Gauss map of the discrete approximate BKL dynamics through successive values of 
u [16]. During an velocity-term-dominated phase of the evolution between "collisions 
with curvature walls," the real part of the speciality index is approximately constant 
and the imaginary part is approximately zero, and the time interval between such 
collisions corresponds to a single point in the discrete dynamics. One can thus directly 
compare the 4-dimensional information in the numerically evaluated speciality index 
with the sequence of discrete BKL values labeled by the number of iterations n, as 
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Figure 8. A plot of the speciality index for the transient true spike solution of 
Lim for w = 1.5 and a; = 10 typical of the case x ^ away from the location of 
the spike at a; = 0, with a closeup showing the formation of the second pulse near 
the origin for «; > 1 which grows relative to the first pulse as w increases. In the 
direction towards the singularity, the first pulse begins at the right most Kasner 
point, making a large clockwise circuit of the origin passing through the left most 
Kasner point, followed by a tiny counterclockwise circuit of the origin coming to 
rest at the intermediate Kasner point. 



well as determine the Taub time interval between the collisions. Since both speciality 
indices are four-dimensional, we can establish a direct correspondence which allows one 
to easily estimate some of the relevant parameters (the p^'s and the relation n ^ At) 
from iS from the numerical evolution, giving an alternative procedure with respect to 
the one implemented in the past by Berger [6, 7]. 

5. Gravitational waves? 

To some extent, interpreting the Blanch! type IX spacetimc as a closed gravitational 
wave, no matter how convincing the arguments by Wheeler and Grishchuk, 
Doroshkevich and Yudin and reviewed by King [4], who fills in the details omitted 
by the others, is a debate more about semantics than substance. The common lore in 
general relativity is that gravitational radiation is associated with the Weyl curvature 
tensor, but when no asymptotically flat background is available for comparison (or no 
short wavelength perturbation of a cosmological model is being considered), it is not 
really clear how to characterize what behavior in the Weyl tensor can be interpreted 
as gravitational radiation. All of these discussions rely heavily on special forms for 
the spatial metric in analogy with the two independent polarizations of gravitational 
waves on less ambiguous background spacetimes. 

However, the speciality index in the context of Mixmaster dynamics seems to be 
an ideal candidate for an object that is independent of all 3-dimcnsional quantities 
except the choice of time in terms of which it is expressed and which can be argued 
in some sense tracks a gravitational wave evolving in time. In the asymptotic BKL 
dynamics, a single collision with a curvature wall generates a complex pulse in this 
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quantity in between intervals of constant real values which change with each collision, 
as modeled by the exact Bianchi type II vacuum Taub solution discussed in Appendix 
B. In the generic Mixmaster dynamics far from the initial (or final) singularity where 
the universe point is simultaneously interacting with all the curvature walls, each one 
is generating overlapping such pulses which form the profile of the speciality index 
graph versus time, but in the approach towards the initial singularity, one has a series 
of temporally isolated pulse transitions between approximately constant behavior as 
illustrated in Fig. 5, like the pulses in an electrocardiogram. 



6. Concluding remarks 

After an entire book [2] has been devoted to the spatially homogeneous Mixmaster 
dynamics including the views of many experts in advanced numerical and qualitative 
techniques, and even more sophisticated analyses have followed in the intervening 
years, one might have wondered what more could be said on the subject. However, 
bringing the speciality index into the discussion, originally introduced in the widely 
different context of modern numerical relativity applied to gravitational radiation 
problems, adds an interesting new perspective to a long and familiar story, one which 
continues even at present as numerical [29] and theoretical techniques advance for more 
general inhomogeneous spacetimes where the speciality indc;x is yet to be exploited 
[17, 20, 23]. Although our focus here is primarily on the spatially homogeneous 
Mixmaster universe, the exact inhomogeneous Gowdy example found by Lim shows 
that the speciality index gives a gauge-invariant lens through which to view the spike 
properties of this solution, and as such has promise as another important tool for 
understanding more complicated inhomogeneous numerically-generated spacetimes. 



Appendix A. Bianchi type IX geometry 

The spatially homogeneous orthogonal dual frame 1-forms in terms of which the 
metric is simply expressed have the local coordinate form (useful for computer algebra 
calculations) ui^ = dt and 

uj^ = cos d^ + sin ip sin ^ d(^ , = sin tpdO — cos ip sin 6d(j), = dV' + cos 6 dcf) 
and the orthogonal frame itself as eo = dt and 

ei = — cot 6 sin + cos ^89+ ^^^^ dj, , 

sm 

62 = cot 6* cos 3^ - sin V' 90 + £^£^ 5 , ^ 63 = , 

sm^ 

where the local coordinates have the ranges = t G [0, 00), x^ = tjj G [0,47r], 
a;2 = 6* G [0, tt] and = (j) G [0, tt]. 

Rescaling the spatial members of the frame and dual frame yields the natural 
orthonormal frame {e^} and dual frame {co"} 

[60,61,62,63] = [9t,a"^ei,6"^62,c"^63] , [a;°,a;\a;^,a;^] = [dt,au}'^,bw'^,cu)^] . (A.l) 

This in turn has a natural associated complex null transverse [30] Newman-Penrose 
frame [14] 

' " 71^^° ' 71^^° ~ ' = "^(^2 + « 63) . (A.2) 
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Letting a dot denote differentiation with respect to t, the vacuum Einstein field 
equations for this metric in Ricci form reduce to 



(dbc) ' 
abc 

{abc)' 
abc 



1 



0, 



{abc)' 
abc 



1 



2\2i 



:r-{c'-a') 



0, 



2a^b^c^ abc 



a consequence of which is the so-called Hamiltonian constraint 

(In a) -(In b) ' + (In b) '(In c) " + (In c) '(In a) ' - ^ (a* + 6^ + c^) + ^ (a^b^ + ?,2c2 + c^a^) = . 

The nonvanishing spin coefficients in the natural Newman-Penrose null frame are 
given by 



1 



P= -M= - 



2^/2 



c b .a 
c b be 



A = -C7 = - 



1 



2^2 



c_b 
c ~ b 



2V2 



a I 
a 2 



cab 
ab be ac 



(A.3) 



and the nonvanishing Weyl scalars are 



^0 = ^4 = ^ 



a2-c2 



ab ac 
ab ac a^b"^ 



2abc 



3c)-2^(a 



V2 = i 



ac 2 ab 
ac be ab 



.a 



4^62c2' 



(A.4) 



a262 



a^h^c^ 



Aabe 



-2aa + ^(a2 + ^2 - c?) + -(a2 - ^2 + c2) 
c 



The Weyl scalars enter the Bianchi identities, which in this case (with ■0o = V'4> 
etc.) imply Stj^Q — 5*1^0 — ^"4^2 — 3*1^2 — 0, and 

Dtp2 = -AV'o + 3pV2 = AV'2 , DtpQ = {p- 4e)V'o - 3AV'2 = A^-o • (A.5) 

Introducing the scale invariant ratio z = V'2/V'o, one finds from the previous Bianchi 
identities 



Dz = Az, Dz = X{3z^ -l) + 2{p + 2e)z. 
Explicitly the second of these is 



1 



abc 



i3z'-l) 



, be . 
In ^ +1 



2a' 



abc 



(A.6) 



(A.7) 



Of course all of these equations are easily re-expressed in terms of the logarithmic 
metric variables and the Taub time. 

Finally the two scalar invariants of the Weyl tensor whose ratio defines the 
speciality index have the following expressions in terms of the Newman-Penrose 
curvature quantities 

/ = VoV'4 - 4^)1 V'S + 3^)2 , J = V'oV'2V'4 - V'l ^"4 - V'oV'3 + 2V'lV'2V'3 - ^"2 • (^-8) 
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Appendix B. The Taub Bianchi type II vacuum solution and Bianchi type 
IX bounces 

The metric of the exact Bianchi type II vacuum solution found by Taub [31] in the 
Taub time gauge is given by Eq. (13.55) and Table 8.2 (p. 107) of [14], namely our 
Eqs. (2.1), (2.2) with 

(w^w^a;'') = {dx^ - x^ dx^ , dx^ , dx^ ) , diu^ = uj'^ A uj'\ dw^ = = dw'' , 

{a,b,c) ={X-\Xe^^,Xe'^^), (B.l) 

where 

= fc-icoshfcr, k = 2{AB)^/'^, (B.2) 

and the initial cosmological singularity abc — > occurs at r — > — oo. By defining 
e = sgn A = sgn B and 

K = (AB)i/2 = k/2 , «oo = e(B/A)i/2 ^ ^ 3) 

one has 

/p2KT , -2A:t\ 1/2 

A = K/u^ , B = Kuoo , X = i ^ j . (B.4) 

Note that simple transformation Uoo has the effect of interchanging the second 

and third directions, which is a symmetry of the spatially homogeneous frame and 
hence of the dynamics. The transformation Uoo —Uoo is equivalent to the time 
reversal transformation r — r. The value Uoo = 1 corresponds to the locally 
rotationally symmetric Petrov type D case, where the speciality index is identically 1. 
The value Woo = is clearly not allowed. 

With these new parameters one has the Kasner limits 

[a, 6, c] [(4if )i/2eT^^ (4if)-i/2e^(±i+"~ )^ (4i^)-i/2eif(±i+»»)r] (3 5^ 

from which one can read off the asymptotic Kasner indices from the coefficients of r 
in the exponentials divided by their sum 

, - _ iS'(=Fl,±l+U~\±l + Moc.) _ (- + 1, U±oo{u±oo + 1)) 

(Pi,P2,P3)±oo = ^77— r- — rr- ^ = 7- ■. — 2 • 

-ft:(±l + lioo + Uoo) 1 + Wioo + Uioo 

Thus comparing with Eq. (4.4), u±oo are seen to be the asymptotic values of the 
Lifshitz-Khalatikov parameter u for this solution, which represents a single bounce 
against a curvature wall W\ in the Hamiltonian picture as illustrated in Fig. 6 of [9], 
under which this parameter simply changes sign. The locally rotationally symmetric 
type D Kasner solution corresponding to Uoo = corresponds to asymptotic motion 
parallel to the curvature wall, which is not possible. 
The metric variables themselves are then 

[a,6,c] [{AKf/'^eP'^±\ {4.K)-^'^eP-'^^\ {'iXy^'^eP'^^^^] , 

ahcdr (4if)-V2e*±Mr = e*±^-^"«^^)'''*±)d(5±T) = e^±df± , (B.6) 

where 

5±=± K{1 + + u^^) , f± = 5±T - \n{{4.KY/H±) (B.7) 

is an afRne transformation from the Taub time to the canonical Taub time of the 
corresponding Kasner limits at infinity. 
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When < A < B then u^c G (IjOo) hes in the interval in which the Kasner 
indices arc ordered: pi < P2 < Pz- This is tlic initial state for the bounce in 
the reversed time direction towards the initial singularity. To reorder the Kasner 
indices after the bounce one then does an interchange of the axes which maps 
u_oo — — "oo —(1 + u-oo) = Uoo — 1 as described in [9]. The asymptotic value 
of the speciality index then undergoes the real shift 



l-Ur,c,\ / 1 + M„o + M 



1 + Woo y V 1 - Woo + w^ 



(B.8) 



on the Kasner interval S e [0,1], which is a big jump in the speciality index on 
this interval when Uoc is just slightly greater than 1, illustrated in Fig. 1. The value 
u = 1 corresponds to the locally rotationally symmetric case of normal incidence to 
the curvature wall. When u^o increases to the golden mean value w = 1/(m — 1) = 
(l + \/5)/2 w 1.618034 which is a fixed point of the BKL map [32] . the self-intersection 
of the circuit descends to the horizontal axis as shown in Fig. 2, and then disappears 
for larger u values. When u reaches the value 2, the right horizontal intercept moves 
to the right end of the Kasner interval [0,1] as shown in Fig. 3. A typical situation for 
u > 2 is shown in Fig. 4. Large values of u^o lead to small changes in the asymptotic 
values of the speciality index according to Eq. (B.8). 

Following the evolution of the speciality index through the bounce one sees 
that the real part roughly undergoes a single large deviation away from the pair 
of asymptotic values, while the imaginary part undergoes a double peaked S-shaped 
departure from the axis, with its midpoint zero roughly correlated with the extreme 
value of the real part pulse. This pair very roughly looks like a single extrema away 
from constant asymptotic values on either side and its derivative graph. The particular 
values of the two independent parameters (X, Uqc) of the Taub solution compress or 
stretch these pulse profiles horizontally (in time) or vertically (in value) but leave their 
correlated shapes roughly qualitatively the same. In fact the complex speciality index 
curve and hence the extreme values of its real and imaginary parts and magnitude are 
completely determined by the Kasner parameter Uoc, while the curve's parametrization 
in time is determined by the bounce parameter K, which stretches/compresses the time 
parametrization, thus changing the location in time and width of the pulse profiles in 
time. 

The explicit expressions for the individual Weyl scalar invariants / and J are 
manageable but not very enlightening. One finds 

I = 4^e-6(»oo+«i )^Vi(woo, C,S), J= __e-4(«oc+«^ )^V2(woo, C, S) , 



2 f3{Uoo,C,S) 



^ = < 7}: r7y (B.9) 

/4(Woo,0,i) 

where fi are complex polynomials and C = cosh(2i^T), S = sinh(2_fi'T). Thus the 
parameter K only rescales the Taub time parametrization of the speciality index 
curve in the complex plane, so this curve (i.e, its image) only depends on the Kasner 
parameter Uqo- However, the plots of the real and imaginary parts of the speciality 
index versus r, although they all share the same shape corresponding to K = 1, are 
stretched or compressed horizontally by this rescaling parameter. In particular the 
locations in time of the pulses and their widths in these plots are stretched/compressed 
by this parameter compared to the value K = 1. 
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The Taub Bianchi type II speciality index satisfies the following symmetries 

^{u^,t) = ^{Uoo,t) , S{-Uoo,t) = S{uoo,-t) . (B.IO) 

The first is an invariance under reflection across the normal line to the curvature wall, 
while the second extends to a local symmetry the asymptotic relation 5(— Woo,oo) = 
S{uix, — oo) satisfied by the real limiting Kasner formula under the reflection u —u 
across the curvature wall direction which results from the bounce with the wall. 
Alternatively the combined change (moc, ''') ^ {^Uoo, —t) leaves the metric invariant 
while changing the orientation used in the duality operation, leading to the complex 
conjugate of self-dual objects. 

In the context of Mixmaster dynamics approaching the initial singularity, the 
temporal pulses in the speciality index are isolated in time. Fig. 5 illustrates this 
behavior for the Mixmaster universe corresponding exactly to r e [0, 200] in Fig. 3(a) 
of the Hobill review in Hobill et al [2], during which three wall collisions take place 
in the recoUapse phase of the evolution towards a big crunch. The initial data 
for this numerical solution are (a(0), /3(0), 7(0)) = (0,0,-1), (a'(0),/?'(0),7'(0)) = 
(—1,0.1,0.3417); the numerical solutions were easily obtained using the computer 
algebra system Maple. Each collision results in a circuit of the speciality index in 
the complex plane, revealing itself as a pulse of correlated temporally isolated pulse 
transitions between approximately constant baselines in the real and imaginary parts 
as functions of the time, separated by Kasner phases of the evolution, leading to a 
profile in time resembling an electrocardiogram. 

Appendix C. Lim's true spike Gowdy solutions 

Exactly this paired pulse- like behavior is seen in the expansion-normalized variables in 
recent work on Gowdy solutions by Lim [17] (see section 4.5), except that they occur 
as a function of the spatial inhomogeneity coordinate rather than of the time. In fact 
these pulses occur both in time and space as a manifestation of a tubelike complex 
pulse surface present in the speciality index plotted versus one of the two nontrivial 
coordinates of the Gowdy spacetime. The Gowdy line element is 

ds^ = e(^-3^)/2dT2 - e(^+^)/2da;2 - e^'-^dy + Qdzf - e'^'M^^ , (C.l) 

where r is a sign-reversed Taub time gauge time coordinate for which the initial 
cosmological singularity occurs at r — > 00. 

As an example, consider Lim's exact Gowdy spike solution of the vacuum Einstein 
equations of section 4.5 of [17] (with his Eq. (34) for Q here corrected by dividing its 
right hand side by 4), in which his alternative Kasner parameter w = 2u+ 1 appears 

P = 2t + ln(sech(wT)) — ln[l -|- {we^sech{wT)x)^] 

- ln(20o) , 

Q = - Qow[e"^'' -|-2('u;tanh(wr) - l)a;^] Q2/4 , 
A = - 41n(sech(wT)) + 21n[l + (we'^sech(wr)a;)^] 

- {w^ + 4)t + A2 . (C.2) 

With the parameter choice w = 0.1, Qo = 1, ^2 = 1 and Q2 = 0, this solution belongs 
to the class < ti; < 1 of "permanent true spikes," which exhibit a single curvature 
transition and hence are most similar to the spatially homogeneous Taub solutions of 
Appendix B. The initial (r — > —00) and final (r — > 00) asymptotic Kasner states in the 
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approach towards the initial singularity common to all values of the spatial coordinate 
a; 7^ in the single BLK curvature transition correspond here to the parameter values 
{w + 2, w) = (2.1, 0.1), for which the speciality index has the corresponding asymptotic 
values given by Eq. (4.5) to be {S-^,S^) « (0.7716,0.9704). Instead for a; = 0, the 
initial and final Kasner states correspond to (w+2, 2 — w) ~ (2.1, 1.9) and the speciality 
index makes a real hyperbolic tangent-like transition in time from iS_oo ~ 0.7716 to 
»Soo ~ 0.6369, leading to a straight line segment along the real axis in the complex 
plane connecting these two Kasner points. All the constant t curves in the singular 
speciality index surface pass through this line segment connecting the two components 
for which x > and x < 0. This discontinuity in the speciality index surface at x = 
caused by the discontinuous limit of 5 as r ^ oo towards the initial singularity is 
the "true permanent spike" lending its name to this class of solutions of the Einstein 
equations, and its location is characterized by the vanishing of the imaginary part of 
the speciality index for all time. Fig. 6 illustrates this situation. 

One can plot the complex speciality index horizontally versus either the vertical 
a; or r coordinate to obtain a tubelike surface with a gap between the two vertical 
lines representing the common asymptotic Kasner values of the initial and final states 
after the inhomogeneous curvature wall bounce. Fig. 7(a) shows a plot of the one 
bounce speciality circuit in the complex plane (horizontal cross-sections) as a function 
of the spatial inhomogeneity coordinate x along the vertical axis. Each circuit of the 
speciality index from the common asymptotic initial and final Kasner values (the two 
vertical lines through the real axis) corresponds to a temporally isolated correlated 
pulse pair in the graph of the real and imaginary parts versus time, in contrast with the 
pulses revealed in the expansion- normalized variable pair (A^_, Sx) graphed versus x 
(see below). Fig. 7(b) shows a plot of the speciality index in the complex plane 
(horizontal cross-sections) instead as a function of the Taub time coordinate r along 
the vertical axis for a sequence of constant values of x. Each curve represents a 
constant value of the time coordinate, which undergoes a complex pulse during a short 
time period which changes with the value of x. In each case the real and imaginary 
parts of the tubelike surface plot reveal the pulse behavior, versus the time or spatial 
coordinate. 

Note that plotting X versus Y for a pair of real functions of a single variable is 
equivalent to plotting their complex combination X + iY in the complex plane. This 
combination A''_ +iT,x of spatial connection and extrinsic curvature in the expansion- 
normalized metric variables (defined by Eq. (13) of [17]) is reminiscent of the Ashtekar 
variables [33], and their individual profiles in space shown at the bottom of Lim's Fig. 7 
[17] or even in Figs. 1-4 of Rendall and Weaver [18] reveal the same isolated pulse 
pair behavior as the temporal profiles of the real and imaginary parts of the speciality 
index, thus revealing a circuit in the complex combination iV_ -|-iSx- Similar circuits 
occur in the general discussion of Uggla et al [23] , in plots of some of the (expansion- 
normalized) real quadratic Weyl scalar invariants against each other shown in their 
Fig. 9. However, the speciality index plot has the advantage of forcing the Kasner 
phases of the evolution to sit on the compact interval [0, 1] of the horizontal axis. 

The remaining Lim "transient true spike" solutions w > 1 exhibit two such 
curvature transitions, and hence the speciality surface has two pulse components, 
distinguished from the permanent true spikes by having common initial and final 
Kasner limits for all x. For both permanent and transient spike solutions, the 
speciality index is completely insensitive to the so called "frame transitions" exhibited 
in the expansion-normalized variables, which are a manifestation of the choice of 
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spatial coordinates. Fig. 8 shows a typical speciality plot for constant x for the case 
w — 1.5 which reveals the small second pulse which forms as w increases past the 
value w = 1 and grows relative to the first pulse as w continues to increase. The 
speciality index during the initial (r — > — oo), intermediate and final (r ^ oo) Kasner 
phases corresponding to the values {w + 2,w,2 — w) = (3.5, 1.5,0.5) has the values 
{S-oo,Sint,Soo) « (0.9635,0.2915,0.4424) calculated from Eq. (4.5). Details wih be 
studied elsewhere. 
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